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Abstract 

We calculate expressions for the state-dependent quasiparticle lifetime, the thermal conductiv¬ 
ity K, the shear viscosity g, and discuss the spin diffusion coefficient D for Fermi-liquid films in 
two dimensions. The expressions are valid for low temperatures and arbitrary polarization. In 
two dimensions, as in three dimensions, the integrals over the transition rates factor into energy 
and angular parts. However, the angular integrations contain a weak divergence. This problem 
is addressed using the method of Miyake and Mullin. The low-temperature expressions for the 
transport coefficients are essentially exact. We find that ~ TlnT, and ~ for arbitrary 
polarizations 0 < R < 1. These results are in agreement with earlier zero-polarization results of 
Fu and Ebner, but are in contrast with the discontinuous change in temperature dependence from 
InT at R = 0 to at 0 < R < 1 that was found by Miyake and Mullin for D. We note that the 
shear viscosity requires a unique analysis. We utilize previously determined values for the density 
and polarization dependent Landau parameters to calculate the transition probabilities in the low¬ 
est order = 0 approximation”, and thus we obtain predictions for the density, temperature and 
polarization dependence of the thermal conductivity, shear viscosity, and spin diffusion coefficient 
for thin ^He films. Results are shown for second layer ^He films on graphite, and thin ^He-^He 
superfluid mixtures. The density dependence is discussed in detail. For k and g we find roughly 
an order of magnitude increase in magnitude from zero to full polarization. For D a simialr large 
increase is predicted from zero polarization to the polarization where D is a maximum (~ 0.74). 
We discuss the applicability of ^He thin films to the question of the existence of a universal lower 
bound for the ratio of the shear viscosity to the entropy density. 

PACS numbers: 67.30.E-, 67.30.ep, 67.30.hr 
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I. 


INTRODUCTION 


Fermi-liquid theory, developed by Landaitii^ in the mid-1950’s, showed how low-temperature 
collective excitations and thermodynamic properties of strongly-interacting normal many- 
fermion systems could be encoded in a few parameters, the Landau parameters, and that 
these parameters were related to a certain limiting value of the microscopic scattering func¬ 
tion.- In Ref. y, Landau also introduced a kinetic equation to describe the nonequilibrium 
properties of a Fermi liquid. The kinetic equation is of the same form as the classical Boltz¬ 
mann equation with the local quasiparticle energy ep(j(r, t) playing the role of a Hamiltonian. 
The application of the linearized Landau kinetic equation to the calculation of transport 
coefficients for bulk ^He has been very successful. In this manuscript we shall apply this 
approach to a strongly interacting many-fermion system in two dimensions. Reviews of 
the bulk calculations at zero polarization can be found in the works of Abrikosov and 
Khalatnikov,- Pines and Nozieres,- and Baym and Pethick.- The calculation of transport 
coefficients for Fermi liquids in three dimensions with arbitrary polarization can be found 
in Anderson, Pethick and Quader-, and Meyer ovich.- The former set of authors used a slick 
general notation that emphasized the similarities in the calculations of the various coeffi¬ 
cients. There exist some measurements of transport coefficients as a function of polarization 
for bulk ^He. Bun, Forbes, Puech, and Wolf^, and also Akimoto, Xia, Adams, Candela, 
Mullin, and Sullivart^ studied the shear viscosity. Sawkey, Puech, and Wolf— studied the 
thermal conductivity. 

Abrikosov and Khalatnikov (AK) in particular showed that the integrals involved in the 
collision integral factor neatly into a product of integrals involving angular variables and 
those involving energy variables. The resulting expression for the kinetic equation could 
then be brought into the form of a linear integral eigenvalue problem for essentially the non- 
equilibrium part of the fermion distribntion function. The exact solutions of these integral 
eigenvalue problems are derived, and discussed in detail by Sykes and Brooker— and also 
Jensen, Smith, and Wilkins.— 

In recent work, we have utilized the kinetic equation approach to study the transition 
between collective excitations in the ballistic regime (zero sound) and collective excitations 
in the hydrodynamic regime (hrst sound) in thin, arbitrarily polarized Fermi-liquid films.— 
For sound, the kinetic equation is usually solved by rewriting the integral equation as an 


2 


(infinite) set of algebraic eqnations by using a Fourier expansion, and then taking moments 
with respect to the angular functions. This procedure is not unique, and we have compared 
and discussed in detail the predicted propagation speeds and attenuation for two different 
approaches.— In the above cited works we have utilized previously calculated^'I^ density 
and polarization dependent Landau parameters in order to obtain numerical predictions for 
thermodynamic and collective excitations for the specihc case of ^He dims. In this paper we 
shall use these same Landau parameters to calculate predicted values for the density and 
polarization dependent transport coefhcients in thin ^He dims. 

The calculation of transport coefhcients for thin ^He dims has been considered by Fu and 
Ebner,— and also by Miyake and Mullin.— Fu and Ebner applied the variational approach 
that was developed by Baym and Ebner— in order to calculate transport coefhcients for 
^He in superduid ^He bulk solutions. The variational approach of Fu and Ebner as applied 
in two dimensions does not lend itself to analytic solution, nevertheless, they were able to 
extract the lowest order temperature dependencies together with numerically determined 
coefhcients for the thermal conductivity k, the drst (or shear) viscosity rj, and the spin 
didusion coefdcient D all at zero polarization. Fu and Ebner obtained In T behavior for two 
of the coefhcients, and pointed out that the source was a weak divergence in the momentum 
space integrals. 

Miyake and Mullin (MM) derived an exact expression for the spin didusion coefdcient for 
two-dimensional fermions with arbitrary polarization. They indicated that in two dimensions 
one obtains a logarithmic divergence at finite temperature in one of the angular integrals if 
one proceeds by strictly following the three-dimensional approach developed by AK. They 
identided the source of the divergence at dnite temperature as an artidce of using zero- 
temperature values for the Fermi momenta in an integrand of one of the angular integrals in 
the kinetic equation. In a very clever analysis, by generalizing the analysis to low but dnite 
temperature they were able to extract an expression that yielded a logarithmic divergence 
only in the zero-temperature limit. In Sec. [TTl we shall derive this fundamental result in 
detail. 

In Sec. iniwe apply the MM method to calculate the state-dependent quasiparticle lifetime 
at arbitrary polarization. This calculation is similar to that of the transport coefhcients but 
is simpler. This allows us to utilize the MM approach in a clear context. We shall compare 
the present result for the quasiparticle lifetime to a previous onel^ that was obtained using 
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a method that is completely independent of MM. In Secs. IIII Al and IIII Bl we calculate the 
thermal conductivity, and the shear viscosity, respectively. In Sec. IIII Cl we include only a 
brief summary of the calculation of the spin diffusion coefficient since that transport coeffi¬ 
cient was analyzed in detail by MM. We note that as in three dimensions the calculations of 
the thermal conductivity and the spin diffusion calculation are very similar. However, unlike 
three dimension, for two dimensions we find that the analysis for the shear viscosity needs 
significant modification. As in the case addressed by MM, the problem in the shear viscosity 
calculation is identihed as being due to the incorrect use of the zero-temperature limit in 
the integrands of the angular integrals. In Sec. IIVI we utilize Landau parameters that were 
previously determined for second layer ^He films on a graphite substrate, and also for thin 
film ^He-^He mixtures to compute density, temperature and polarization dependencies for 
the transport coefficients. Our results for the shear viscosity are used to calculate the ratio 
of the shear viscosity to the entropy density. Sec. |V] is the conclusion. 


II. QUASIPARTICLE LIFETIME 

We examine a system of N = + A|, spin-i fermions in a box of area L^. The particles 

have bare mass m, and interact with two-body potential V (r) that is assumed to depend only 
on the scalar distance between the particles. The particles hll two Fermi seas up to Fermi 
momenta and k^, and we introduce the convention that the spin-down Fermi sea will 
always be the minority Fermi sea in the case of nonzero polarization. The term polarization 
denotes the magnetization per particle which will be denoted by V, thus V = M/N = 
{N-f — Ni) /N. The terms coverage and areal density {N/Li^) are used interchangeably. The 
system is assumed to be at some finite but low temperature T in the sense that T « Tfp 

The quasiparticle lifetime due to quasiparticle-quasiparticle interactions in two-dimensional 
Fermi-liquids was calculated in Ref. Q. The method used in that reference was borrowed 
from two-dimensional electron theory, and took advantage of the similarity in structure 
between the collision integral and the free fermion dynamic structure function. The 
fluctuation-dissipation theorem together with Stern’s analytic expression^ for the two- 
dimensional susceptibility yielded an analytic expression for the low-temperature lifetime. 
In this section we shall repeat this calculation using the Miyake-Mullin method. This is 
convenient because the lifetime calculation is similar but simpler than that for the transport 
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coefficients. The results from this section are in agreement with our previous results, and 
will be used in the following sections for the transport coefficients. 

After some simplification,- the quasiparticle collision frequency is given by: 

rh-EEE fT(l, 2, 3, 4)(5(ei T £2 *^3 ^4)<5p^+p2,p3+p4^(Ti+0'2,cr3+cr4^2^3^4 > (^’l) 

^ P2,0'2 P3,o-3 P4 ,o-4 

where ripo- = l/[exp(/3(epo- — p ■ u — jUa-)) + 1] is the Fermi distribution function, /3 = l/ksT, 
ficr is the chemical potential for the cxth Fermi sea, and u is the fluid velocity. In this section 
we can set u = 0. The hF’s are transition rates, and we have defined 


n 


pa 


— 1 flpa 


1 4- e-^(ep<T-p u-pa) 


( 2 . 2 ) 


The standard treatment in three-dimensions follows Abrikosov and Khalatnikov,— and in¬ 
troduces new integration variables in terms of energies and angles. These integrations are 
independent of one another, and in lowest order in temperature one can find a closed form 
expression for 1/r in terms of an angular average of the transition rates. Label the incoming 
quasiparticles as Pi,P 2 and the outgoing quasiparticles P 3 ,P 4 . The standard angular vari¬ 
ables, 6 and cp, are defined as follows: 6 is the angle of p 2 measured relative to the direction 
of pi, and (f is the angle between the planes formed by the pairs of vectors {pi,P 2 } and 
{p 3 ,P 4 }. As discussed by MM, in two dimensions <4 can only take on two values: 0, tt. We 
illustrate these two possibilities in Figs. [T] and |2l Along with MM we shall refer to these 
two processes as forward and backward scattering, respectively. We note that the forward 
and backward scattering processes have a direct and exchange relationship since Fig. |2] is 
obtained from Fig. [T]by exchanging ps and p 4 . 

Figs. [U and [2] define the conventions that we shall use throughout this paper to label 
the angles associated with the quasiparticle momenta. All angles are measured counter¬ 
clockwise. The angles with i = 1,2, 3,4 are the angles of pj as measured from pi -|- p 2 . 
The angle 6 ij is the angle of p* as measured from the direction of pj. In the discussion 
below, we shall find it convenient to use the following definitions: = On, a = 6 ^ 43 , and as 

noted above 9 = ^ 2 - 

Since the interaction is assumed spin-independent, total spin is conserved in the collisions 
as indicated by the Kronecker delta in (12.Ij) . Thus, we must have <73 = ai and a 4 ^ = a 2 - 
In the spin parallel case (T 2 = cri the exchange diagram is identical to the direct diagram, 
and therefore they must be counted only once in the phase space integrations. In the spin 
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FIG. 1 . The momentum space diagram for the forward scattering process, tt < 9 ^ < 2 tt. From 
momentum conservation pi + P2 = Ps + P4- The angle 0 is the angle of p2 as measured from the 
direction of pi; 6 ^ is the angle of ps as measured from the direction of pi + P2. We shall also 
need: the angles <^3 = ^31 and <^4 = @41 which are the angles of ps and p4 as measured from the 
direction of pi, respectively; a = ^43 is the angle of p4 measured relative to ps. 


anti-parallel case ai = —a2 the direct and exchange diagrams give different contributions to 
the total transition probability. Performing the spins sums, replacing the momentum sums 
by integrations, and performing an integration over p4 yields: 

^ j dp2dp3 + ^2 - 63 - £ 4 )^ 2^3714 , 

1 1 

-1- 

^(Ticri ^(Tl—(T1 

where for later use we have introduced spin parallel and spin anti-parallel collision frequen- 
cies. In this expression we have set A = 1 . Thus, the units of the are energy-time“^- 

area^. As usual, the factor of one-half appearing with the spin-parallel transition probability 
prevents over counting as discussed above.- 


( 2 . 3 ) 
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FIG. 2. The momentum space diagram for the backwards scattering process, 0 < 0$ < tt. We 
note that this figure can be obtained from Fig. [1] by exchanging p 3 and p 4 . From momentum 
conservation pi + P 2 = Ps + P 4 - 


We first consider the spin-parallel lifetime, and separate ont the angular integrals: 

^ ^ + ^2 - es - e4)n2n3n4 , 


Tr 


CTiai 


2 ^oo roo 

J P2dp2 J Psdps 


- 27 r 


dO 


dOs Wp (9)6 (ei + 62 — es — e4)n2n3n4 , 


( 2 . 4 ) 


where we have taken advantage of the symmetry in 6 about vr, 63 is defined as the angle of 
P3 measured with respect to Pi + P2, see Fig. [H and the subscript f or b on W identifies the 
transition probability as that for forward scattering (ti < 63 < 27 r) or backward scattering 
(0 < 03 < tt), respectively. We now rewrite the angular variable 63 in a more useful form. 
With an eye on Fig. [T]the law of cosines yields: 


pl=pl + f - 2p3i cos 03 , 
003 = 


dpi 


2 p 3 £sin 03 ’ 

where following MM the useful variable i has been dehned: 


( 2 . 5 ) 

(2.6) 


^ = |Pl + P2 


( 2 . 7 ) 
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It is now convenient to introduce the angle a = 6^43 as shown in Fig. [T] 

— sin 63 = PzPa sin a . ( 2 . 8 ) 

From the law of cosines again: 

=pI+ pi - 2 p 3 P 4 cos (tt - a), (2.9) 

p\+pi + 2 pip 2 cos 9 = pI+ pi + 2 p 3 P 4 cos a. (2.10) 

Eq. fl2.10p can be simplihed using energy conservation: 

£1 + £2 = ^3 + £4 , (2.11) 

where we have dehned Si = pf/2m*. The quasiparticle label on the effective masses is needed 
since at hnite polarization the effective masses are state dependent. We End: 


P 3 P 4 cos a = {ml - ml) ^ 3 + piP 2 cos 6 , 


( 2 . 12 ) 


where we have dehned the important quantity ^ 3 : 

c - pI-pI 

2 ml ’ 

where ^3 ~ OlksT). We use fl2.12l) to eliminate a in 


( 2 , 13 ) 


plpisiv? a = pIpI sin^ 9 + (mlpl — m^pf — PiP 2 {'ml — ml) cos9) ^3 — {ml + m’D'^^l. (2.14) 


Finally, we combine Eqs. fl2.6p . fl2.8p . and (12.141) to yield: 
(i6*3 =- 


where 


dpi 


2p3P2xlsm‘^ 9 - —^3 - I ^3 

ei2 pfpi 


(vr < 6*3 < 27r), 


(2.15) 


1 _ A ^ P2 A I A ^ Pi .A 1 

— = 1 H-cos 9 ] -1 H-cos 9 — 

ei2 \ Pi J e2 \ P2 / ei 

.. PF2 1 A , Pfi a 

1 H-cos 9 -1 H-cos 9 

Pfi J (■F2 V PF2 J ^fi 


(2.16) 


The F subscripts on the Fermi energies and Fermi momenta indicate that we only need the 
zero-temperature limit for 642 since .^3 itself is 0 (/cbT). 

















We note that Eq. fl2.15p is exact. For the spin parallel case the linear term in ^3 vanishes 
making the thermal correction term ~ 0 ((/cbF)^). The thermal correction changes from 
linear order in ^3 for anti-parallel spin scattering to quadratic order in ^3 for parallel spin 
scattering. With Eq. fl2.15p . Eq. fl2.4p becomes: 


' 0 - 1(71 


1 

¥ 


dp2 dp3 dpi 


dO 


sin2 e-i^ 


Wp + €2 — €3 — ei)n 2 n^ni . (2.17) 


The integrals are brought to their hnal form by introducing dimensionless variables Xi = 
- p)- 

1 




dx2dxzdx^6{xi + X2 — X3 — X4)n2n3n4 


dO 




sinH 


It is convenient to split the 6 integral into pieces: 


t-A 


de = 


dO + 


dO + 


de 


t-A 


( 2 , 18 ) 


( 2 , 19 ) 


where following MM we have dehned A = l/(/9p). At low temperature A 1 and so the 9 
integral can be approximated as 


r-^ de 

sin 6 * 


+ 


T —A 


de 


02 


\2 
£3. \ 


( 2 . 20 ) 


The integrand in the second term is not singular, and so the second term can be neglected 
relative to the hrst since it does not contribute to lowest order with a logarithmic temperature 
dependence. Thus we obtain the MM form for the angular part of the two-dimensional 
collision integral: 

1 2{mlJpkBTf 


r—A 




h^P. 


dx 2 dx 3 dxi 6 {xi + X 2 — X 3 — X4)n2n3U4 


( 2,21 

sm e 


The energy integrals are evaluated by Morel and Nozieres^ and the hnal result is 7 r 2 / 4 . 
The transition rates can be Fourier analyzed as usuaU^ yielding: 




( 2 . 22 ) 


e=Q 


where the T^posO) = cos {£9) are Chebyshev polynomials of the hrst kind,— and the pa¬ 
rameters (Uq = 1 and ai = 2 for i > 1. We can now introduce the lowest order ‘T = 0” 
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approximation by replacing the full transition rate by its £ = 0 value ~ . The 

remaining angular integral can now be evaluated: 


T —A 


d9- 


sin 9 


2 In 


2ei 

ksT 


The hnal result is: 


1 


Taiai 2 Cl 


261 

keT 


(2.23) 


(2.24) 


The calculation for anti-parallel spins proceeds analogously. From Eq. fl2.3p : 

2 /*00 /*00 f*7V 






dp 2 


dp3 ‘^^(5(ei 62 - 63 - e4)u2n3U4 , 




> p7T . p'K 

d€2 6^63 / d9y d9^ 


r-27r 


d93 Wp 


(2.25) 


(2.26) 


where — is for forward scattering vr < ^3 < 27r, and -|- is for backward scattering 0 < ^3 < vr. 

\2 roo 


X 5(61 62 - 63 - 64)n2U3n4 . 

We now use fl2.15p for anti-parallel spins: 

dds = 


±dpl 


2piP2x klV? 9 - ;^^3 


1 2 m*(m*y 


hViP-fTi Jo 

pTT 

X / [ Wp -^^ + Wp -'^^] 


d62d63(i64(5(6i -|- 62 — 63 — 64)?7.2U3n4 

d9 


sin^0-^e3 


(2.27) 


Performing the energy integrals yields: 

1 PmlPmpV 


T — A 




{ksTy 


1 + wp-^^] - 


d9 


sin 9 


(2.28) 


where in this case A = pksT /| 6 i 2 |. Utilizing the £ = 0 approximation we obtain: 


(i.3r)2,,T^) , (2,29) 


2^4 Pa^P-. 


0-1 


ykeT ) 


By combining fl2.24p and fl2.29p we obtain the total quasiparticle-quasiparticle collision fre¬ 
quency at hnite polarization: 


TT 


{mypWppHn 


261 

ksT 


[m 


Pa, 1 


P-a, 2 


- {W’y’' + In 


4|6i2| 


{ksTf 


6l 

(2.30) 

In this expression the momenta and energies are zero-temperature Fermi momenta and Fermi 
energies (see also Eq. fl2.16p h 
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1. Zero polarization 


At zero polarization, Eq. fl2.24p for the spin-parallel collision frequency is still valid. 
However for the spin anti-parallel collision frequency l/ei 2 = 0 in the zero polarization 
limit, and thus the leading order correction in the denominator of dO'i is quadratic. For zero 
polarization Eq. fl2.29p becomes: 


1 

2 h^ 




a —<j 

f,0 







(2.31) 


By combining this result with fl2.24p we obtain the total zero-polarization collision frequency: 

26 p 


i=l(p = 0) 

To Ta 




2 h^ 




keT 


(2.32) 


where the tq notation will be used below. The same quantity in Ref. [ij Eq. (3.24) differs 
by the appearance of a factor of 3/8 in the leading coefficient instead of 1/2. This is not 
a problem since in that derivation the coefficient of the log term is uncertain with regards 
factors of 0(1) because of the vagaries of the low-temperature limiting process. We note 
that one cannot obtain fl2.32p by taking the zero-polarization limit of fl2.30p . 


S. Full polarization 


In the full polarization limit one simply omits the contribution from the anti-parallel 
spins in fl2.3Up : 


1 

R 


i(p = i) 

R 


2h^ 6^t 


2ei7’-|-\ 

1^) ■ 


(2.33) 


III. TRANSPORT 


The derivation of the transport coefficients in a two-dimensional Fermi liquid proceeds 
in a very similar way to that in three dimensions. Thus, this and the following sections on 
transport coefficients will necessarily be brief. For the details we refer the reader to Baym 
and Pethick- for example. We shall concentrate on those aspects that are specihc to hnite 
polarization and two dimensions. The general transport equation can be written: 




T p^pcr 


Ve 


per 


rper , 


(3.1) 
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where hpo-(r) is the local quasiparticle distribution function dehned with the local quasi¬ 
particle energies epo-(r) . As usual the local quasiparticle distribution function is expanded 
around local equilibrium-^ 

^pcr(r) Up(7 [ep(j(r)] -|- (5?7,pCT(r), (3.2) 

where is the Fermi distribution function. Equivalently, we can expand the local quasi¬ 
particle energies around a set of local equilibrium energies epCT(r): 

epa(r) = ^p<7 + Y1 • (3-3) 

p' 

The collision integral on the right hand side of the transport equation can be written; 

^ ^ hF(l) 2,3,4)(5pj_|_p2jp3+p4(5(j]^+o-2,o-3+o'4*^(^i T ^2 ^ 4 ) 

P2.P3.p4 (3.4) 

X [hih2(l - h3)(l - 714 ) - (1 - hi)(l - h2)h3h4] . 

As discussed in Sec. m the sums over p 3 and p 4 include only distinguishable final states. We 
now expand the collision integral fl3.4p to first order in the 6np^{r). Consider the products 
of distribution functions in the square brackets of fl3.4p . and substitute fl3.2p . This yields: 

[•••] = -13 (Cl + C 2 - Cs - C 4 ) nin 2 n 3 n 4 , (3.5) 

where we have defined 

dvi' 

= 0^“^* = -(3niniCi , (3.6) 

0€i 

and we have made use of the identity: 


[?7,i?7,2n3n4 - nin2n3n4] (5(ei €2 - €3 - q) = 0 . 


(3.7) 


Performing the spin sums the collision integral becomes: 

P2,P3,P4 ^ ^ 

(Cl + C2 - C3 - Cl) ninirijni, 


(3.8) 


where to lowest order we have now replaced all of the remaining local energies e by local 
equilibrium energies e. As discussed in Sec. HIl the value of 6*3 determines whether the 
transition rate is given by the forward or backward scattering diagrams in Figs. [U and 121 
respectively. Eq. fl3.8p is the starting point for the collision integral for each of the transport 
coefficients discussed below. 
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A. Thermal conductivity 


The thermal current is given by 

jth — 'y ^ 

a 

where 


dk 

7^ 




Suka = TT^Vk^ ■ V(/csT)g^(eko 


de 


kcr 


(3.9) 


(3,10) 


and Eq. fl3.10p defines the function q^j. With no loss of generality, we can assume that 
'^{ksT) is in the x-direction. The y-component of jth vanishes, and fl3.9p becomes: 


X r°° 

{jth)x = j^J dkkY^ 


dn 


k(T 


r‘27r 


de 


kcr 


(eica - {kBT)\q„{e]^^)v 


Fa 


d(p COS^ 0 , 


h2 


T|V(T)|5^ 


* 2 

^a'^Fa 


Oti 

dx-^q„{x)x. 


(3,11) 


By definition the thermal conductivity is given by jth = —kVT, thus we obtain: 

dn 


K = 


■ 


T 


dx^ q^{x) X. 
dx 


(3.12) 


The integral over q^x can be evaluated exactly from the kinetic equation by using the 
results of Sykes and Brooker.— In Eq. fl3.8p we have expanded the collision integral fl3.4p to 
linear order in the nonequilibrium distribution functions 5npo,(r). Now we need to proceed 
to do the same for the left hand side of (13.11) . Since the system is assumed to be in steady 


state we can immediately set 


dh 


pa 


dt 


= 0. We systematically ignore the gradients of 6 n- 


po-( 


Then with u = 0 and no polarization gradient we have: 

dn 


Vh^a ^ ?^Ve 


de 


pfT 


p, + ^vr(r), 


— ~ _ 
^p^pcr — 


de 


p^pcr ? 


(3.13) 

(3.14) 


pa 


and thus the left hand side of 


To lowest order VpCpo- ~ 


C/T) 

becomes: -^|^VT(r) • Vp^ 
^pCpo- = Vpcr is the Fermi velocity. 


= /9(ei-/ii)^(VT)-vi. 


We note that we have omitted a V/i contribution in the expansion. This term, propor¬ 
tional to the entropy, was shown by Sykes and Brooker to contribute only to the even part 
of q, and that this is of higher order than the odd part. Below we shall point out that our 
go-(a:)’s have odd symmetry. 
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Using fl3.6p and fl3.10p we find: 


Ci = Vj ■ V{kBT)q^{ei). 


(3.15) 


The kinetic equation can then be written: 


/3(ei-yUi)nini(vi ■ VT) = ^ ^pi+P2,P3+P4<^(ei + ^2 - £3 - e4)ni?7,2n3n4 

P2,P3,P4,o-^ 


x(V">"'.(9) [(vi . Vr)g„(ei) + (v 2 ■ Vr),,;(e 2 ) - (vs ■ VT)q„(e3) - (v4 ■ Vr)g,;(e4)] . 

(3.16) 


where it is understood that in the case of spin parallel scattering the range of the 63 in¬ 
tegration is restricted from 0 to vr. It is straightforward to show that the angle between 
Vi and VT is arbitrary. Define 7 such that vi ■ VT = ni|VT|cos 7 . Then similarly 
for the other quasiparticle velocities: Vj ■ VT = Vi\'VT\cos9i = nj| VT| cos ( 7 -|- 6 *ji) = 
ni|VT| (cos 7 cos 6^11 — sin7sin0ji), where i = 2,3,4. As a reminder, from Figs. [U and [ 2 ] we 
have 621 = ^2 = 9, 6*31 = $ 3 , and 9^1 = $ 4 . We note for further reference that we will need 
the following expressions in our analysis of backward scattering at 0 < T < 1 : 


cos $3 
cos $4 


2p2 sin^ 9 


cos 9 + 


£2 

2pip2 sin^ 9 
£2 


sin $3 
sin $4 


sin 9 
sin 9 


2 p 2 (pi +P 2 cos 9) 

jpl - pI) 

£2 


(3.17) 


where € = pi -|- p 2 as dehned in fl2.7p . The integral fl3.16p is symmetric in 9 with respect 
to the intervals [0,7r] and [tt, 27r]. Thus the terms containing sin 0*1 must cancel out. From 
(I3.17P it is clear that this is true for any value of 6 * 3 . The term |VT| cosy cancels out from 
both sides. We then divide through by Vi to obtain: 


{^1 ~ Pi) niUi — ^pi-i-p2,p3+p4<5(ei -|-€2 — €3 — e4)?T.in2n3n4 


P2,P3,P4,cri 




V2 


9(71 (ei) H-cos 6*219^' (€2) 

Vi 1 

- cos 93 iq^^ (€3) - — cos 94 iq„f (64) 
Vl ^ 


(3.18) 


We now transform the momentum integrations into energy and angular integrals as was 
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done for the collision time in the previous section: 


2 /‘°° 

(3 (ei - /ii) nirli = J de2de^dei5{ei + €2 - €3 - e4)nin2n3rr4 


Pai J 0 




sin^ 6 — ^ 


<?< 7 i (ei) + ( 62 ) cos 6 > - (es) - (€ 4 ) cos 9 


w * )2 .. pp-i-i 


(^) 


P(TlP—<Jl Jo 

,* '\2 /-TT 


Sin 


20 _ 


U-, 


U-n 


Qai (ei)H- —q-ai (^ 2 ) cos 9-q^^ (€ 3 )-(€ 4 ) cos 6 ^ 




P(TiP—ai Jo 


sin^ 6 ^ — (—' 

' £12 ' 


V _. 


V-r 


In (ei)H- —q-<,i (£2) cos e-g„ (€3) cos $3 - —q-„ (€4) cos $4 


(3,19) 


where we have used the symmetry in the integrals about 9 = tt, and ^3 and ei 2 are dehned 
in Eqs. fl2.13p and fl2.16p . respectively. At this point we can show that to lowest order in 
temperature the terms in fl3.19p that are proportional to cos0 can be omitted. For either 
parallel or anti-parallel spins we have: 



Wl^,'i9)cos9 

^ sin 2 0-(^)2 


^/Ao 


d9^^ 

sm 9 


0 . 


(3.20) 


For the spin parallel case: A = ksT/ep^ <C 1, and for the spin anti-parallel case: A = 
■\//csT/ei 2 1. We note that from (12.161) the quantity €42 depends explicitly on 9\ 


1 

^12 


1 

€2 



2 (m2 — ml) 
P1P2 


cos 9. 


(3.21) 


This expression would seem to introduce a problem at P = 1. However at full polarization 
these terms are removed by the minority-spin Fermi velocities in fl3.19p . and make no con¬ 
tribution to the physics in that limit (see, for example Fq. (I2.33p ). At hnite polarizations 
0 < P < 1 the 9 dependence in 642 introduces different values of A at the lower and upper lim¬ 
its. Nevertheless since these differences are of 0(1) or less we still have A = ^/ksTje^ 1 
at both limits. 

Following the same line of argument as in three dimensions, we introduce dimensionless 
variables: Xi = jJ^ei — /ij). If we then let Xi —Xi, we now see that qa{xi) is an odd 
function of its argument: qa{xi) = —qa{—Xi), in lowest order of temperature for the thermal 
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conductivity. If we let —X 3 and X 4 —)■ —X 4 , we obtain: 


ximrii = —{ksTY 


X 


(m 


o —00 
* \3 /•TT 


jp. 

t'(T\ 


I dO — , i-^^— / aa — , ' 

0 ^sin^ 6 — P<TiP-ai Jo ^sin^ 9 — (^) 


dx2dxodxii 5{xi + X 2 + x^ + x^)nin2n3n4 

wp^^io) 


<A’n>_„r fpg_ wp-«m 




* \2 fTT 


/ de 

PaiP—ai Jo 


PuiP—ai Jo 

wp-'^pe) 


d9 




Qaiixi) + Qaiixs) 




sin^ e 
P 


daiiXo)- 


V-„^ sin^ 9 


V. 


Cl 


P 


(l-ar{x4) 


(3.22) 

At this point we can relabel the variables X 3 ,X 4 X 2 due to symmetry under the integral. 
Further, we note from Eq. (I2.30j) that the first set of angular integrals is just the collision 
time: 


4 1 2 


P T, 


{ml) 


* 'vS /-TT 


d9 


WJp9) 


^2 f-TT 


Pa Jo 


sin2 0- (^)2 

WrP9) 


.. + <^ld9 


WJ-‘^{9) 


PaP-a Jo 


sm 


20 _ (i^' 

^ ei 2 ^ 


d9 


PaP-a Jo 


siP9 


?3 


(3.23) 


We follow MM (see Sec. IIII Cl below) and introduce generalized frequencies: 


4 2 

—z/, = j-{kBT) 


TT" V P~ 




2^ UmP)^ wrpo) 2pl^siiP9 


PaP-a Jo 


siP9-(^] 


^ PaP-a Jo 


sm 


^9-(^ 


P 


V-a 2p^p_a- sin^ 9 

Vrr P 


(3.24a) 


. (3.24b) 


For simplicity in notation we have switched from ai to a. Then the kinetic equation becomes: 


/ OD 

dx 2 dx 3 dx 46 {xi + X 2 + Xo + X4)nin2n^n4 

■00 


4 1 


4 / n-o 


— — qa{xi) + qa{X 2 ) - :pl^aqa{X 2 ) + ^ U-aq-a{X 2 ) \ ■ (3.25) 

TT^ To- L J TT^ Tl^ \ n„ J J 

In the third term in curly brackets we have introduced the notation for the areal density of 

the Fermi sea, and this should not be confused with the similar looking Fermi distribution 

function (12. 2 h whose subscript is a momentum label and not a Fermi sea label, thus: 


n„ = 


A ■ 


(3.26) 
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The energy integrals can be found in Appendix A of Bayni and Pethick:- 


dx^, dxi 6 {xi + X2 + Xy, + X4) ^ 377,4 = 


Xi + X 2 
X — Q-{xi+X2) ’ 


f°° j j , r/ \ 1 X 4 + TT^ 

/ dx2 dx3 dXid[xi + X 2 + x^ + x^j n2n3n4 = -- — . 

J-00 2 1 + e 

Thus, the kinetic equation can be brought into non-diagonal Sykes-Brooker form: 


(3.27) 

(3.28) 


71 

Y I = 


dX2K{xi,X2) qcxixi) - (1 - VaTo) qa{,X2) - ] q-a{x2) 


n-n 


rin 


(3.29) 


where the Sykes-Brooker kernel is defined as: 

^ ^ (1 + e-^ 2 )(e^ 2 -a:i _ 1 ) 


(3.30) 


This kinetic equation mixes the two components of q^- It is in very similar form as the 
kinetic equation for spin diffusion as derived by MM. In Sec. IIII Cl we shall briefly write 
down the relevant expressions for the sake of comparison. We rewrite fl3.29p with a matrix 
representation of the coefficients of qa{x 2 ). 


/ OO 

dX2K{xi,X2) qajxi) - y^^Xaa'qa'{x2) 

•OO ^ t 


(3.31) 


The coefficient matrix is given by: 


( 1 — Z/fTf \ 

. (3.32) 

1 -nn/ 

The matrix diagonalization proceeds by using the general method described in Anderson, 
Pethick, and Quader.- We note that in this case A is not symmetric. The eigenvalues of A 
are: 

A+ = 1, A_ = 1 - (z/tTt + . (3.33) 


The ± subscripts on the A’s refer to the plus and minus roots of the quadratic equation 
generated by diagonalizing fl3.32p . In a spinor sense, plus and minus label the top and 
bottom rotated spin state, respectively. 

We introduce transformed variables f = St and q = Sq where the transformation matrix 
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S and its inverse are given by: 


5 = 


% n-t 

n'Tt + nn 


(3.34) 


= 


n^Tt 

nti^t 


We find: 


r = 


i^fTt + ly^T^ 
1 


q = 


^tn(^ + 

fi4, fit 

ng+ _|_ ng- 
% % 


\ ^tg+ _ ffjg- 
\ H % 

In terms of these variables the diagonalized pair of kinetic equations are: 


vr _ 


dX2K{xi,X2) qa{xi) - \aqa{x2) 


(3.35) 

(3.36) 

(3.37) 


As shown in Eq. fl3.12p . the important quantity is not itself but rather the integrated 
quantities: 

roo 

(3.38) 


dn _ 

Qa = -J dx—qa{x)x. 


From Sykes and Brooker Eq. (60) the solutions to the diagonalized problem can be written: 


Qa = 


TT^fk 




(3.39) 


2(3-A.)' 

where H{\) is an infinite series involving the eigenvalues \± that will be explicitly written 
below. We substitute fl3.36l) into fl3.38l) and then inverse-transform using: 

Qt\ _ ln^^iQ+ + niT^Q. 

Qi ] - n^nQ- 

The exact solution for the low-temperature thermal conductivity in two dimensions can be 
written: 

vrfc? 


(3.40) 


K = 






T ^ , 


where 

Qt = 


71 


Z/fTt + 2 L 

1 71^ 

Z/fTt + 2 


nn(i + KM)) 

I^trt(l + (ntM)) 


MK 

3-Ah 

H{K 

3-Ah 


+ (z/|Tt - (n^/nt)z/;r^) 

+ inn - in/H)nn) 


HjX-) 
3-A_ 
HjX.) 
3-A_ 


(3.41) 

, (3.42a) 

n, (3.42b) 
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(3.43) 


together with the explicit expression for H{X): 




4n + 5 


4 ^(n + l)(2n + 3)l(n + l)(2n + 3)-A| 


According to Eqs. (I2.30p and (13.24^ we hnd i^o-To- ~ 1/lnT. Thus, at very low temperatures 
we can set A_ ~ A+ = 1. In this limit then Eqs. fl3.42p simplify to Q„ = ^H(l)Ta , and the 
thermal conductivity becomes: 


K, = 




Trr 


(3.44) 


1. Zero polarization 


At zero polarization the eigenvalues are A+ = 1 and A_ = 1 — 2z/ro. By inspection of 


fl3.42p we have Qf = Qi = (7r^/4)if(l)ro and therefore: 


TT 


niV ^ 0 ) ^ ^kgTepH(l)Ti) , 


(3.46) 


in agreement with (I3.44p . In the £ = 0 approximation this becomes: 


k{V = 0) = vrh^/c 


B- 


Hil) 


{eF/ksT) 


m* 


-I- IE' -I- IE 
f,0 ^ f,0 ^ '^'^6,0 


rt I 


tt 


In {2eF/kBT) 


(3.46) 


2. Full polarization 

At full polarization all quasiparticles are in the t Fermi sea, and thus the terms with 
the spin anti-parallel contributions z/^i- and do not appear. By inspection of the kinetic 
equation fl3.3ip the eigenvalue A = 1. From fl3.42p we have = (7r^iJ(A)/2(3 — A))ri and 
therefore: 

f{V = 1) = *:|r£„ff(l)ri, (3,47) 

in agreement with fl3.44p . 

Summary 

We hnd that the temperature dependence for the thermal conductivity at arbitrary po¬ 
larization 0 < P < 1 is K~^ ~ TlnT. This is in agreement with the zero-polarization results 
of Fu and Ebner.— 
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B. Shear viscosity 


We consider a Fermi-liquid film flowing with speed Ux in the x-direction. The flow is 
not uniform. There exists a small non-zero ^/-gradient of the velocity dux/dy that will drive 
a transverse momentum flux axy The coefficient of proportionality rj is the hrst or shear 
viscosity: 


(^xy y' 


dUx 

dy 


(3.48) 


The stress tensor axy can be written in terms of the non-equilibrium part of the distribution 
function:- 


^xy — 


E 


dp 

7? 


Px 


de 


pa 


dpy 


6 n 


pa 5 


(3.49) 


we note that Vpa^y = dep^/Opy. Then using fl3.6p we can write 5npa in terms of the driving 
field: 

fin ■ fin ■ 1 fi'n 

(3.50) 


^_dni driil dux , . 

OTli — 2 {Pixl^iy T Piy'^ix) Qa y 


where from symmetry: PixViy = PiyVix = m*VixViy. For the shear viscosity we will need the 
hrst two terms of (I3.50p in powers of (epo- — p). Following Sykes and Brooker we obtain: 


^xy 


[E/^fe 


dn 


x„ 


x:^pa,y} 


^'"(1 -f^)mlv^^xVa,yqaiepa) 

per P^F<J 


1 dux 


dy 


(3.51) 


where = (3{epa — p) was introduced before Eq. (I2.18p . The shear viscosity is then: 

'PipFaVpaY [ ^(1 + ^^)cos^ 6 >sin^ 6 '|^g<^(ep^). (3.52) 

^ J jSepa oepa 


P = 


We now perform the angular integration, and change integration variable from p to x\ 


TT 


P = 


4h2 


E 


(PFa)"^ Xrr . dn 


m: 


dx (1 -f- 


[depa dx 


qa{x) 


(3.53) 


We can resolve qa{x) into symmetric and antisymmetric components: qa{x) = qY\x) -f 
qf\x). Because of the even and odd symmetry of the two components of the integrand, we 
can write p = -|- 7 ^^“^ where: 



P 


(a) 


TT 

IF 



dx^qfix) 


TT {PFa)"^ 1 

4h2 m* (depa 



dx x^q^J^\x) 


(3.54a) 

(3.54b) 
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As in the previous section, the integrals over can be evaluated exactly.— It is straight¬ 
forward to show that ~ O(T^), the same as in three dimensions, and thus we can 

ignore the contributions of the anti-symmetric part of in the remainder of this discussion. 


In the absence of thermal or polarization gradients, the left hand side of the kinetic 

^^pfj \ 1 r T dUa 


equation fl3.ll) in leading order reduces to:— 
fl3.50p . the kinetic equation becomes: 


de 


per 


2 [PxVpa,y + PyVpa,x] 


Then with 


fl\'fl\V\yP\x ^ ^ ^Pi+P2iP3+P4^(*^l A (^3 64)?T.i?7.2?7.3?7.4 

P2,P3,P4 

X [ViyPi^qaiiei) -f V2yP2xqai{e2) - (^s) “ (^ 4 )] 

+ wp~'"^{e) [ViyPi^q^^{ei) + V2yP2xq-aA^2) - V3yP3xqai{(^3) ” V4yPixq-aA(^^)] 

+ [ViyPi^q^^{ei) V2yP2xq-aA(^2) - VsyPSxqaii^s) - V4yPi^q_^^{ei)] I , 

(3.55) 


where we have canceled out a common factor of dux/dy. This is similar to the kinetic 
equation for the thermal conductivity fl3.16p except that in this system both the x and y 
directions play special roles. Thus we need to include information as to the directions of 
the momenta with respect to the x-direction, say. We introduce angle 7 which is the angle 
between pi and the x-axis: pi • :£ = pi cosy. 

At this point we need to emphasize a key difference between the thermal conductivity 
calculation, and that of the shear viscosity. For the shear viscosity we need to permit the 
incoming and outgoing quasiparticle momenta to differ slightly from the zero-temperature 
Pf^s in order to obtain sensible results. Of course, energy and momentum still must be 
conserved in quasiparticle collisions. This situation is also discussed by Novikov^i in his 
treatment of the shear viscosity for a two-dimensional fermion system. 

In principal all four momenta can be unequal to the zero-temperature ppi^i = 1,2, 3,4. 
In our model however we shall fix pi = ppi and p 2 = Pf 2 i and only permit ps and p 4 to differ 
from their zero-temperature values. This simplifies the calculation while still maintaining 
the essential features of the the effect of hnite temperature. In addition, this maintains 
consistency with the MM treatment of the divergence in the integral over 9 in the kinetic 
equation as discussed in Sec. [TTl In Fig. [3]we illustrate the sort of scattering process that can 
yield a nonzero value from the phase space integral that appears in the kinetic equation. In 
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FIG. 3. Momentum space diagram illustrating the effects of letting the quasiparticle momenta 
move slightly off of the zero-temperature value pp for the backwards scattering case. In this zero- 
polarization example we set pi = P 2 = Pf as discussed in the text. It is clear that $3 / 9, and 
also since pi and p 4 are not parallel $4 7 ^ 0 . 


terms of the angular variables, Eq. fl 3 . 55 p becomes: 


niTii sin 7 cos 7 = E <^P 1 -|-P 2 ,P 3 -|-P 4'^(^1 + ^2 ~ ^3 ~ 

P2,P3,P4 

X I wp { 9 ) sin 7 cos (ei) sin (7 9 ) cos (7 9 )qai (€2) 

- sin (7 + cos (7 -h ^pq^Pes) - sin (7 -h <h{) cos (7 <h{)gai {ep 
+ Wp ( 9 ) sin 7 cos 'jq^^ (ei) sin (7 9 ) cos (7 (62) 

- sin (7 -h cos (7 -h $3)9^1 (63) - sin (7 -h cos (7 -h <h{)g-ai (q) 


wp 


{ 9 ) sin 7 cos 7g^i (ei) -h sin (7 9 ) cos (7 9 )q-a^ (62) 

sin (7 $3) cos (7 -h ^3)qai (es) - sin (7 $4) cos (7 -h ^’Pq-a^ pp 


( 3 . 56 ) 


where for convenience we have dehned = {v-aP-P/paPa). 

Expanding the trigonometric functions, we obtain terms like sin 7 cos 7 cos( 26 *) and 
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sin0 COS0 cos( 27). The terms proportional to sin 6 ' give no contribution to the summa¬ 
tion due to invariance under the transformation 6 —)■ —6. The right hand side of the 
equation can be cast into a product of energy and angular integrals as was done in the 
previous section, and 03.561) becomes: 


^2 2 


niUi = (ksT) — / dx2dx^dxi6{xi + X2 — X3 — X4)nin2n3n^ 

• ^ ./ —rv^ 


(m 


* \3 PIT 


X 


Pai Jo 


de 




sm^e 


•_C 3_'|2 - 


QaPxi) + (lai(,X2) COS (26*) 


- Qai (xs) COS (2${) - {x^ cos (2${; 




PcriP—cri Jo 


sin^d- (^) L 


(xi) -h (xa) cos {26) 


wp-^pe) 


PaiP-ai J0 ^silPO - (^) 


- Qai {X 3 ) COS (2${) - (a; 4 ) cos (2<h{) 

(xi) (xa) cos {29) 

■ (xs) cos ( 2 $^) - (X 4 ) cos ( 2 <h^) 


(3.57) 


The second and third integrals in Eq. 03.571) describe forward and backward scattering 
between spin anti-parallel quasiparticles, and their exact forms depend on polarization. Ap¬ 
propriate expressions will be derived below where we will discuss three polarization ranges 
separately. At this point we can evaluate the contribution from spin-parallel forward scat¬ 
tering. The calculation of cos<h{ and cos<h{ begins by determining cos 6^3 and cosd 4 . From 
Eq. ([23]); 


Pi + P 2 p + Pi 

2 P 3 IP 1 + P 2 I 
Pi + P 2 I ^ pI-pI 

2p3 2p3|pi-hp2 


(3.58) 


We note that the factor of |pi -|- P 2 I in the denominator of 03.58P is the source of the 
enhancement of the role of head-on scattering as will be discussed below. 
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We shall begin with the spin-parallel channel. Thus as per the above discussion we can 
set Pi = p 2 = Pf- From fl2.13p we have the following dehnition: 

& = (3.59) 

Using conservation of energy 2p|, = p‘i+ p‘i we identify: 

pI-pI = • (3.60) 

Thus since |pi - 1 - P 2 I = 2pi?cos (|), we can write to lowest order in T: 

cos(^ 3 ) = cos(|) + g , (3.61) 

2 cos(|) 

where we have dehned A 3 = Then exchanging ps and p 4 in fl3.58p we have: 

cos( 04 ) = cos (f) - g . (3.62) 

2 cos(|) 

Eqs. (13.6ip and (13.621) can now be used to determine $3 and $ 4 . In the case of forward 
scattering between parallel-spin quasiparticles, $3 and ${ can be signihcantly different from 
0 and 9 respectively only at 0 ~ tt, head-on scattering. Thus using fl3.6ip and fl3.62p . one 
has 


sin<F{ = sin( 6'3 - 6 »i), 


and similarly 


Therefore 


and 


sin 9i cos 6*1 — sin 9i cos 63 
tan (|) A 3 


. ./ . . tan (f)A 3 

sm $4 sm 6 *-^- 


cos(2<F3) 1 — 


tan^ (I) A 


tan^ (I) A 


cos( 2 ${) cos(26') — 


where we have used 6*1 = 27r — Substitute Eqs. (I3.65p and 


fl3.57p . and one obtains: 


[m 


* \ 3 


r-A 


p2 

Fai 


dO 




(^) 


sin 9 




(3.63) 

(3.64) 

(3.65) 

(3.66) 

into the hrst integral in 


1 — tan^ ( - 1 A 


(laAx2) 


(3.67) 


Note we have assumed the parity of go-(x) and the change of variables have already been 
applied. 
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1. P = 0 Zero Polarization 

At zero polarization, we set pf\ = Pfi, and = 1- The contribution from the 

forward scattering of spin-antiparallel quasiparicles is identical to that of the forward scat¬ 
tering between spin-parallel quasiparticles as analyzed in fl3.67p . Thus, the second integral 
in fl3.57p is identical to fl3.67p with the substitution Wj"" . 

For the backward scattering between spin-antiparallel quasiparticles, one notices that 
03 ^ 271 — 01 and 04 ^ 0i as 0 is not close to vr, therefore 


and 


Then 


and 


sin $3 = sin( 6'3 — 0i), 


sin 01 cos 01 — sin 0 i cos 0 ^ , 
tan (I) A 3 


= sin 0 -|- 


sin $4 = sin( 6*4 — 0i) , 


sin 01 cos 01 — sin 0 i cos 04 , 
tan (I) A 3 


cos( 2<F3) = 1-2 sin^ $3 , 


. s tan^ (I) A2 
cos(26')-- 


(3.68) 


(3.69) 


(3.70) 


cos( 2$4) = 1 - 2sin^(<F4), 

tan^ (f) A 


(3,71) 


Substituting these back into the backward scattering integral, we obtain for fl3.57p : 

*'(3 P +00 

dX2dX3dX46 (Xi + X2 -\- X3 -\- X4) 711+12713+14 


Tii+ii = {ksT) 

I>TT-A 

/ d0 


2 2 im*y 


Pf J- 

Wp "1 (0) + wp (0) + wp (0) 

sin 6 * 


QaPxi) - 


1 — tan"^ 


A' 


Qai {X 2 ) 


(3.72) 
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As before the equation can be cast into the form: 


TT 


^+oo 


-Tomni = 


dx2dx‘idxi5 (xi + a;2 + xs + X 4 ) nin 2 n^n 4 


q{xi) - (1 - ^o'To 


where tq is given by (I2.32p . and we have dehned a generalized frequency: 




r-A 


wp{e) + wj\e) + w'^\d) 




tt/ 


dO 


sin 6 


tan^ - Ao. 


(3.73) 


(3.74) 


Following the same steps as for the thermal conductivity, Eq. fl3.73p can be cast into Sykes- 
Brooker form: 


71 


^+00 


dX2K(Xi,X2) - A(j(l2)] 


(3,76) 


The solution is then 


Q = - 




2i/o ’ 


with the eigenvalue: 


A = 1 - uqTo , 


and from Ref. 


12 


E 


(4n + 3) 


c(A) ^ 

(1 - A) 4 ^ (n + l)(2n + l)[(n + l)(2n + 1) - A] 


(3.76) 


(3.77) 


If we keep only the zeroth order components of the transition rates, and simply set A 3 = A 
to extract the correct order of temperature dependence in z/q-, we obtain 


<^0 « {Wjl + Wji, + Wlj, 

The expression for the zero-polarization viscosity follows from fl3.54p : 

7(P = 0) = A 


(3.78) 


m* 1 — A 

Using 1 — A = UqT ~ 0(1/InT), and c(A 1) ?» |, we find: 

Zh? 

r]{V = 0) = - 


(3.79) 


4 tt irU + Puji ^ 


(3.80) 


The zero-polarization temperature dependence rj ^ ~ is the same as that found by Fu 


and Ebner.— There is an important feature of this expression that needs to be pointed out. 
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Because 1 — A in fl3.79p is proportional to tq the explicit dependence on the quasiparticle life¬ 
time has canceled out. This is quite different from three dimensions where rj is proportional 
to r. 

At this juncture it is convenient to explain why the above unorthodox treatment of the 
kinetic equation is necessary for the fermion shear viscosity in two dimensions. If we had 
proceeded in the usual manner by hxing the Fermi momenta equal to their zero-temperature 
values then we would find no uqTo term in fl3.76p . and thus the eigenvalue A = 1. Unfortu¬ 
nately A = 1 is an eigenvalue of the associated homogeneous integral equation. Since the 
left hand side of fl3.75p is not zero, then there would be no solution. 


0 < F < 1 


In this polarization range, we note that the second term in Eq. fl3.58p is always negli¬ 
gible compared to the hrst term. Thus for the forward scattering between spin-antiparallel 
quasiparticles we set $ 3=0 and = 9. Then the second term in the curly bracket of 
Eq. (I3.57P becomes 


ParP-ai 


de 


wp "'"1 


( 0 ) 


sin2 9 -^ 

ei 2 




(3,81) 


The integral for the spin-antiparallel backward scattering has to be evaluated using Eqs. 
(3.18): 




d9 




( 0 ) 


sin^ 9 -^ 

ei 2 


Qai {xi)- (l-2 sin^ $ 3 j {x 2 )+D^p(^ cos^ 9-cos‘^ $ 4 j( 0 : 2 ) 


(3.82) 

The integral over the last term is much smaller than that for the first two terms. Indeed if 
one sets Wp~‘^^{9) ^ , using Eqs. (3.18) one obtains: 


d9 


cos^ 9 — cos^ $4 


sin2 9 -^ 

ei 2 


= 0 . 


(3.83) 


Thus we can ignore the last term in the integral, and therefore the kinetic equation at finite 
polarization becomes spin-decoupled. Combining Eqs. (13.671) . fl3.8ip . fl3.82p and fl3.83p . the 
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kinetic equation becomes: 

2 r+oo 

niUi ={kBT)'^— / dx 2 dx 3 dx 46 {xi + X2 + Xs + X4)nin2n3n4 


h\ 


p. 


(Tl ^ ^ 

^2 /•TT 


sin 6 ^ 


1 — tan^ 




Qn (^ 2 ) 




PcriP-(Ji Jo 


PcriP—ai Jo 


dO 


2q ?3 
E] 


sin" 0 - ^ 

E 12 


9 ^ 71 ( 2 : 1 ) - qaAx2) 


sin" e-^ 


q<T^ (a^i) - (1 “ 2 sin" $3 j {X 2 ) 


(3.84) 


With the aid of the dehnition of the quasiparticle lifetime fl2.30p the kinetic equation be¬ 
comes: 

vr" ^+°° 

T 

qa^{xi)- 1 - (i/^^ + I/®) 


dx2dxodx45 {xi + X2 + xo + X4) nin2non4 

qaii 


(3.85) 


with the dehnitions 


•'i = T^i-ksT) 


IT 




4 ” = (kBT) 


^ pI JA " sin^ 
tt " „.2m*{m*_^f r-^ 


■ tan 


A 


3 5 




de^ -^sin 2 <|., 


' PaP-cT Ja ~ sin 6 ' 
The kinetic equation can now be cast into Sykes-Brooker form: 


TT 


b+oo 


-Trr = 


with 


The solution is 


dX2K{xi,X2) [qaixi) - \aqaix 2 )] , 

XD 

c(A^) 


Qa = 


2[i'l + uf’ 


(3.86) 

(3.87) 

(3.88) 

(3.89) 

(3.90) 


where c(A) is given in Eq. fl3.77l) . In lowest order of hh, we can obtain approximate analytic 
expressions for the generalized frequencies: 


1/7 


{ksT) 


^tt" {mlf 


p, 


2 '^W.o ’ 


(3.91) 


A3) 


tl dp Q- 


6,0 


Po 


p2 _ p2 

1 + TV- 

‘JPaP-a 


Pa + P-o 


Pa - P-a 


(3.92) 



































We note they are both on the order of T^. 

For the polarization range 0 < P < 1 the shear viscosity becomes: 

_ Svr {PFaY 1 


mz 




where c(A) ~ 3/4 since f^t„ ~ 0(1/InT). In lowest order we find: 
3 


OIL \ "V 






'^F-0 




1 + 


2po-P- 


In 


Pa+P-o 


Per -P-a 


^Fa 


(3.93) 


(3.94) 


Thus, at finite polarization we find rj ^ ~ T^, the same as at zero polarization. Further, the 
quasiparticle lifetime has canceled out in the same manner as at zero polarization. 


3. Full polarization V = 1 


At full polarization we ignore those terms that involve scattering between anti-parallel 
spin quasiparticles since there are no particles in the minority Fermi sea. Thus we ignore the 
^(3) f|3.89p . and therefore A-i- = 1 — f^ti. The solution becomes = c{X^)/{2 f^) 

3 /{ 8 f:^). From Eqs. (I3.93p and (I3.94p the shear viscosity at full polarization becomes: 

dvr (pFf)"^ 1 f 


p{V = 1) = 


32^2 


m. 


t n 


8 vr Wjl yksT 


(3.95) 


Summary 

We hnd that the temperature dependence for the shear viscosity at arbitrary polarization 
0 < P < 1 is ri~^ ~ T^. This is in agreement with the zero-polarization results of Fu and 
Ebner.— 

C. Spin diffusion 

As noted in the Introduction, the longitudinal spin diffusion coefficient D for a two 
dimensional Fermi liquid at arbitrary polarization was calculated by Miyake and Mullin.— 
We shall include a very brief discussion of the calculation of D in order to compare this 
kinetic equation solution to the closely related solutions for the thermal conductivity and 
the shear viscosity as discussed in the previous sections. Further, using our results for the 
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Landau parameters for ^He films we can calculate predicted values for D as a. function of 
density and polarization. These numerical results will be presented in the following section. 

In this case the driving held is a chemical potential gradient. The left hand side of the 
kinetic equation fld.ip can be written: niUiVi ■ (Vpi), and the linearized ansatz for the 
non-equilibrium distribution function becomes: 




(3.96) 


Choosing vi || (V/ii) with no loss of generality, the kinetic equation becomes: 


niUi = E *5pi+P2,P3+P4*^(*^l T 64)?7 .i?7.2?7.3?T.4 

P2,P3iP4 

(6) (ei) -h qai P 2 ) cos 9 - q^^ (ea) - q„^ {e^ cos 9 
+Wp~''p9) gai(ei) + A„^v^^q-„Pe 2 ) cos 9 - q„PeP - A„^q_^PeP cos 9 
+Wp~‘'p9) qaiiei) + A^^q_^Pe 2 ) cos9 - qa^es) cos $3 - A^^q.^PeP cos $4 


(3.97) 


where following the notation of MM we introduce: 


A _ V-ai R-a^/N^ 

— 


—0-1 


•"0-1 


R.JN, 




(3.98) 


The R„ parameters are the proportionality constants that connect the chemical potential 
gradients with magnetization gradients: V/ig- = {R„/NpR))'^ {nV) where n = = 

N/A is the total number density, and the Nq = m*/(27rh^) are the single spin-state density 
of states at the Fermi surface. For completeness, we shall write down the explicit expression 
for the i?o- since this quantity hnds its way into the hnal expression for the diffusion constant: 


Rn = 


crn_„ 


(1 + Fo'-Xl + F^-) 


Tia—a 77’—crcr 




(3.99) 


AT- rr<,(l + Fr + Fo-P/Ns + n_pi + + F^P/N, 

and we associate a = —1} with a = {t, i}, respectively. We note that from the 

requirement of stable equilibrium!^ the numerator of fl3.99p must be positive. The Landau 
parameters are dehned by: 

/•2it jq 

Fr' = NS j —/™; (3.100) 

The reduction of fl3.97p to Sykes-Brooker form follows the same steps as for the thermal 
conductivity. One obtains a kinetic equation in non-diagonal Sykes-Brooker form: 

TT^X /•°° r 1 

— lv = / dx2K{xi,X2)^qaPxi) - {1 - u„^T^Pq^Px2) + i^-aiT„-,q-aPx2)p (3.101) 
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where we have introduced the the spin-parallel and spin anti-parallel scattering times fl2.24p 
and fl2.27p . the generalized frequency 

^2 ^ ^ 


IT 




(3.102) 


and in addition we have also used: 


I'-rr = 


Po 


P-0 


-A„Vrr 


(3.103) 


Eq. fl3.10ip agrees with (B5) in MM which marks the start of their analysis. The only 
differences are in the definition of some parameters: (to-)'^^ = (Tr^/d)^ and {va/fa)^^ = 


VrrTrr. 


We rewrite (I3.10ip with a matrix representation of the coefficients of qa{x 2 )'- 
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To- = 


dX2K{xi,X2) qa{xi) - E Qa' (^ 2 ) 


(3.104) 


where for simplicity we have switched from ai to a. The coefficient matrix is given by: 


A = 


L - -n'B 


(3.105) 


The eigenvalues of A are: 


A+ = 1, A_ = 1 - (i^tn + nn) > 


(3.106) 


where as in Sec. IIII Al the ± subscripts on the A’s denote the rotated up and down spin states. 
The diagonalization is accomplished by transforming to variables f = Sr and q = Sq where: 


= 1 (-n 


(3.107) 


In terms of the transformed variables the diagonalized pair of kinetic equations are: 


0 

T 


’ —00 
/*oo 


dx 2 K(xi,X 2 ) q+(xi) - q+(x 2 ) 
dx2K{xi,X2) q-ixi) - X-q-{x2) 


(3.108) 

(3.109) 


We have written the pair of equations separately in order to emphasize that the equation 
with the unit eigenvalue is indeed homogeneous. 
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As with the other transport coefficients, the important quantity is not itself but rather 
the integrated quantities: 


^ _ r dn^ 
Qa= 


(3.110) 


This definition differs from that of MM by a minus sign. From Sykes and Brooker the 
solutions can be written: 


Q+ 

Q- 


c, 

c(A-) 

2(1-A_)’ 


(3.111) 

(3.112) 


where C is an arbitrary constant to be determined below, and c(A_) is a series involving the 
eigenvalue A_ defined in Eq. (I3.77p . The diffusion coefficient is given in terms of Q, obtained 
from the inverse transform: Q = S~^Q. 

The spin current for the ath component is jo- = \v'^RaQ^)■ MM fixed the value of 
C by assuming that the system is in the frame of reference where the spin current due to 
bulk motion vanishes, thus X^o-Jo- ~ 0- 

The diffusion current is given by j = = —DVm. Thus, -D = | and 

we have: 


D 




( t/tn + nn ^ 


c(A-) 

(1-A_)- 


(3.113) 


This is in agreement with MM’s fundamental result. Using fl3.106p we can also write: 


1 v‘^vlRi;Ri 4n + 3 

2 4 (n + l)(2n + l)[(n + l)(2n + 1) — A_]- 


(3.114) 


The temperature dependence of D depends on whether one is at zero polarization or non-zero 
polarization. These two cases will be discussed in turn. 


1. Zero polarization 


At zero polarization we set Uo- = vp, R^ = —= R and = u in fl3.114p where: 


R 


1 + ^0 
2 


u = 





Wr"{0) sine 

£2 


(3.115) 

(3.116) 
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and A = ksT/ep- The only pole is at the npper limit of the 6'-integral, thus we have 
extended the lower limit to zero. The spin diffusion coefficient reduces to: 


D{V = 0) = 


4(1 + FS) 

8iy 


c(A_). 


In the i = 0 approximation the frequency u is given by: 




*^2TT/o■-(T { ^bT) 


ep 


2,6 p 


(3.117) 


(3.118) 


and the quasiparticle lifetime Tq is defined in fl2.32p . D depends on both the spin parallel 
and the spin anti-parallel transition probabilities through the the eigenvalue A_: 


A_ = 1 — 2vTn = 1 — 2 




ffl ^ f,0 ^ 

Written out in detail, the diffusion coefficient at zero polarization is: 

(1 + 4 ) 1 


(3.119) 


2 VttV 


€F 


X 


(1-A) 


\kBT J 

An + 3 


(3.120) 


(n + l)(2n + l)[(n + l)(2n + 1) - A_] 


This result is in agreement with MM (41). The dependence of D on the spin-anti-parallel 
transition probabilities through the eigenvalue mimics the solution for bulk Fermi liquids.— 
In a more general form than the i = 0 approximation, the eigenvalue A_ depends on the 
transition probabilities hTJand through the angular averages that appear in 

z/o- and To-. 


2. Nonzero polarization 


The low-temperature physics changes qualitatively at any non-zero polarization. The 
important contribution is from the frequency Va (13.102^ : 




2fP-o 


Pa 
2fP-c 


de 




w;^-^(9) sine 


In 


Pa / PaP-a 


P 

Pa + P-o 


(3.121) 


Pa - P-o 


where the integration limits are not cut-off since at hnite polarization there are no poles. 
As pointed out by MM there is a logarithmic singularity in the limit of zero polarization 
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with no singular behavior as a function of temperature. The lifetimes at hnite polarization 
are given by fl2.30p . Thus we find TaVa ~ 0(T“^ln“^(T)). This should be compared with 
fl3.119p above. At low temperature then we have A_ ~ 1 + 0(T“^ln“^(T)). As discussed 
in Refs. 12 and!?, in the limit A —)■ I"*", we can replace the sum by its first term yielding 
c(A) 3/4. 

Thus, the spin diffusion coefficient at P 7 ^ 0 is given by: 

-ikBT)-\ (3.122) 


3 




(I 

2 


m^m’^ 


(Wo) 1° 

Pt+Pi 

Pr-Pi 


The absence of the In (T) divergence in D at finite polarization was hrst noticed by MM 
who pointed that the change in temperature dependence from zero polarization was due to 
the inability of the system to conserve momentum in spin anti-parallel collisions at non-zero 
polarization and low enough temperature. In the limit of full polarization the spin diffusion 
coefficient vanishes. From fl3.122p it can be seen that D{V —?■ 1) 0{pI). 

The finite-polarization spin diffusion coefficient given by (I3.122P does not go smoothly 
into the zero-polarization result (I3.120p . We can easily calculate the value of the polarization 
Vc at which the diffusion coefficients from Eqs. (13.1221) and (I3.120p are equal. In the limit 
of very small polarization |P| <C 1 fl3.122p becomes: 


D{\V\ < 1) ^ - 


3 /hn (l + Fo“) 


8 VttV In \kBT 




(3.123) 


where all terms have been set to their P = 0 values except for the term that is singular 
in that limit. We can set c(A_) = 3/4 in (13.1201) with little error since its range is 3/4 < 
c(A_) <1.— Then the only difference between fl3.120p and fl3.123p are the logarithmic terms. 
By inspection: 

P, = {ksT/AeF ). (3.124) 

For polarizations less than Vc the diffusion coefficient may be measurably larger than the 
zero-polarization diffusion coefficient. Of course this would be within a very small regime 
since this analysis is valid only in the limit ep. 


Summary 

The temperature dependence for the spin diffusion coefficient at zero polarization is 
D~^ ~ In T. This is in agreement with the zero-polarization results of Fu and Ebner.— At 


34 



































finite polarization 0 < P < 1 the temperature dependence for the spin diffusion coefficient 
is D~^ ~ T^. These results were obtained by MM.— 


IV. APPLICATION TO THIN ^He FILMS 


In this section we calculate transport coefficients for thin-hlm ^He systems. The system 
specihc information is provided by the angular integrals of the transition rates [9] that 
appear in the quasiparticle lifetimes r and the generalized frequencies v. The transition rates 
can be written in terms of the scattering amplitudes: {9) = Dimensionless 

scattering amplitudes can be dehned by: 


^(9) = N„a%(e), 


(4.1) 


where Nq = m/(27rh^). The transition rates then become 

Wji {B) =—^\A%[e)\\ (4.2) 

j. m 

In turn, for forward scattering, the Fourier components of the scattering amplitudes can be 
written in terms of the Landau parameters:— 


tt 

«/y = 


tt 

“m = 


//ti+iVoVy) - 

(1 + + Afpy) - ’ 

ftt 

h 


(4.3a) 

(4.3b) 


(1+iVoVy)(i+ivpy) - ' 

We note in passing that in this notation the forward scattering sum rules^^ become hFJ°'(0) = 


0 . 

At zero polarization one can also write the backward scattering transition probability 
in terms of the forward scattering amplitudes:— 


al-^{9) = ay^{9). (4.4) 

It is not known whether a comparable exact result can be obtained for nonzero polarization 
in two dimensions. The important point is that at zero polarization in two dimensions D 
can be written solely in terms of the Landau parameters. This point was made by Miyake 
and Mullin. In fact this is valid for all of the transport coefficients. In order to compute 
backward scattering contributions at nonzero polarization, and also head-on transition rates 
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used for the shear viscosity, we shall proceed by making some reasonable assumptions. For 
backward scattering at nonzero polarization: 




(4.5) 


and for the head-on scattering needed for the shear viscosity: 


= 0 or it ) ss 0^/(0 = it ) 


(4.6) 


In three dimensions one faces a similar problem because one also needs information con¬ 
cerning the 0-dependence of the scattering amplitudes where 0 is the angle between the 
planes formed by the momenta of the incoming and outgoing quasiparticles. An approx¬ 
imate solution at zero polarization was obtained by Dy and Pethick.— Unfortunately the 
s-p approximation does not generalize to nonzero polarization. In two dimensions as noted 
previously 0 can only take on the values 0 and vr. 

For ^He in two dimensions we can calculate Landau parameters to high ordersl^ by 
utilizing effective s-wave and p-wave T-matrix elements determined by experimental data. 
Thus in principle we can also calculate the Fourier sum for WJ^' {9) to high orders. For 
the numerical work to be discussed below, however, we shall use the lowest-order ^ = Q 
approximation for the transition rates. In Fig. 0] we compare Ay and for £ = 0,1 at 
n = 0.0132 A-2 on a graphite substrate. We see, at zero polarization, AJ'*' dominates the 
other three components: as one expects, the singlet channel dominates the 5-wave scattering. 
As the polarization increases however, increases rapidly and eventually becomes the 
dominant component. Therefore we can approximate the transition rates by simply keeping 
the i = 0 components over the whole polarization range. 

We begin by examining the transport coefficients in the forms kT In {2Tp/T), and 
In (2Tp/T) to analyze their density dependence for "P = 0 and V = 1. The results are 
presented in this manner in order to take advantage of the fact that at zero polarization and 
full polarization the explicit temperature dependence factors, see Table [H In Table UTl we show 
the values for the system of second layer ^He on graphite. We include in the table the Fermi 
energies, and the effective masses. The values of the effective masses at full polarization come 


from Ref. 


16 


We can compare the qualitative behavior of the transport coefficients with 
their bulk ^He analogs. Figure 3 in Bedell and Pines^^ shows the pressure dependence at zero 
polarization of kT, and DT"^. In three dimensions each of these quantities appears to 
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FIG. 4. Dimensionless scattering amplitudes (I4.ip versus polarization at n = 0.0132 on a 
graphite substrate. Comparison of with and Ag'*' with A[^. We note that at P = 0 we 
find Ag^ « aI^. For all polarizations AJ"*" > 4A 


be a monotonically decreasing function of pressure. In two dimensions this is not necessarily 
the case. One can easily extract the explicit density and effective mass dependence of the 
transport coefficients by examining Eqs. (I2.32p . (I3.45p . (I3.79p . and fl3.12np . By inspection of 
Table [H we hnd n ~ rj ~ and D ~ Thus, the explicit density 

dependence tries to increase the transport coefficients with increasing density whereas the 
explicit effective mass dependence tries to decrease the transport coefficients with increasing 
effective mass. 

At "P = 0, kT In (2Tp/T) decreases monotonically with increasing density. This behavior 
is dominated by the increase in the effective mass. There is additional density dependence 
carried by the scattering amplitudes. In Table HTTI we include the most important scattering 
amplitudes for the ^He system of Table [TTl The contribution of the scattering amplitudes to 
K is shown in the second column of Table Hm The quantity is non- 
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TABLE I. Expressions for the transport coefficients at P = 0 and V = 1 rewritten in forms that 
are useful for analyzing the density dependence. 


Transport coefficient Polarization Expression 




K 

K 

V 

ri 

D 


V = 0 

V = 1 

V = 0 

V = 1 
P = 0 


Hiiy 


n 


1 


kT In 




|ijf 12 


r]T‘^ = 


n 


1 


V,oi 


Akim‘SJ {m*lm)^ |i^J2 + + |^ttj2 


rjT^ = 


/ Svr/i® 


n 


\k‘yim?J (rn*/m)4 |j;tt |2 
T>r2 In f c(l) 


n 


V.oi 


1 + Ff 


0 


k'^m^J (m*/my 


monotonic, however, the extent of variation is small compared to that of the effective mass. 
This is also the case for the spin diffusion coefficient. The important scattering amplitude 
in this case is which is fairly constant. Thus, the density dependence of the spin 

diffusion coefficient is dominated by that of the effective mass. The case of r]T‘^ is more 
intriguing, as it seems hrst to increase, and then to decrease with density. At low density, 
the cubic density dependence dominates the viscosity. Thus at low densities the viscosity 
increases with increasing n. At higher densities the increase in the effective mass eventually 
takes over, and the viscosity starts to decrease. This “bump” feature is also present in the 
density dependent behavior of Tp. 

In contrast to zero polarization, at full polarization the transport coefficients exhibit 
an increase with increasing density. This behavior can be understood by referring to the 
effective masses at P = 1 shown in Table UTl In the limit of full polarization the effective mass 
shows only a slow increase with increasing density. Thus, at full polarization the explicit 
increases with density dominate the small increases in the effective masses. 

In Tables [TVf IVl and IVTl IVIII we show analogous results for ^He adsorbed to 3.14 A 
and 4.33 A superfluid ^He dims, respectively. It is important to note that these results are 
restricted to a much smaller density range than for ^He on the second layer of graphite. 
The reason for this difference is that in the superfluid ^He environment the ^He undergoes 
a transition to a transverse excited state at an areal density n = 0.036 A“^.— The data in 
Tables lIVtIVIIl cover a density range less than the first three data points in Tables [TTl Hill 
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TABLE II. ^He on a graphite substrate. The zero-polarization V = Q and full polarization V = 1 
thermal conductivity n, shear viscosity 77 , and spin diffusion coefficient D as functions of areal 
density n with the explicit temperature dependencies factored out. The units of kT are (10~® 
erg s“^) ; the units of are (10“® g s~^ mK^); the units of DT^ are (cm^ s“^ mK^). We also 
include the effective masses, the Fermi energies, and we note that D{V = 1) = 0 from its definition. 
The values of k are obtained from (j3.45l) and (j3.47p . for r] from (I3.80p and (|3.95p . and for D from 


(j3.12Up . for zero and full polarization respectively. 


Density (A m* 

/m 

e^ (K) 

Krin(2rp/r) 

rjT^ 

DT"^ ln(2rp/T) 

V = 0 

V = 1 

= 0 V = l 

0 

II 

V = l V = 0 

V = 1 

0 

II 

0.013 1.29 

0.82 

0.52 

1.64 

0.138 

5.4 4.88 

0.38 X 10^ 

5.42 

0.025 1.72 

0.81 

0.75 

3.16 

0.123 

10.4 8.34 

1.41 X 10^ 

2.67 

0.037 2.64 

0.86 

0.72 

4.36 

0.047 

15.2 4.70 

3.01 X 10^ 

0.47 

0.046 3.66 

0.92 

0.64 

5.10 

0.040 

18.7 4.91 

4.61 X 10^ 

0.29 

0.054 4.88 

0.95 

0.57 

5.82 

0.014 

23.3 2.06 

6.79 X 10^ 

0.06 

TABLE III. ^He on a 

graphite substrate. The dimensionless scattering amplitudes g, and 

Landau parameter Fq 

that are 

! the input into calculating the transport coefficients shown in Ta- 

blelE 








Density (A“^) 


+ 1 Vol 

2 + |i 

t4 |2 

6,01 

Fq- 

14'’''*-12 
1^6,ol 

\A^ |2 



0 

II 


0 

II 

0 

II 

V = 1 

0.013 


1.12 



-0.51 

0.64 

0.35 

0.025 


1.44 



-0.62 

0.88 

0.68 

0.037 


1.46 



-0.72 

0.93 

0.79 

0.046 


0.72 



-0.71 

0.46 

0.78 

0.054 


0.90 



-0.79 

0.59 

0.76 


Using Tables IIIIl IVl and I VI II we can compare the Landau parameter Fq, and some of the 
scattering amplitudes for the two substrates. Over the same density range is markedly 
smaller in magnitude in the mixture film than on graphite. The denominators for k and rj 
iVoi^ +1 + iTiP are considerably smaller for the mixture films than for graphite. 
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TABLE IV. ^He in a 3.14 layer film of ^He. The zero-polarization 7^ = 0 and full polarization 
V = 1 thermal conductivity k, shear viscosity ry, and spin diffusion coefficient D as functions of 
areal density re with the explicit temperature dependencies factored out. The units of kT are (10“^ 
erg s“^) ; the units of ryT^ are (10“’^ g s“^ mK^); the units of DT^ are (10^ cm^ s“^ mK^). For 
this mixture film mn = 1.56rre is the hydrodynamic effective mass.— The values of re are obtained 


from (I3.45P and (j3.47p . for ry from (I3.80p and (|3.95p . and for D from (I3.120p . for zero and full 
polarization respectively. 


Density (A 

m* 

P = 0 

Iran 

V = 1 

(-F 

P = 0 

(K) 

P = 1 

reP ln(2rp/T) 

= 0 P = 1 

C-H 

o 

II 

■2 

P = 1 

PT^ ln(2rp/T) 

P = 0 

0.013 

1.31 

0.84 

0.32 

0.99 

0.285 

0.92 

0.098 

0.63 

1.34 

0.016 

1.47 

0.87 

0.36 

1.21 

0.341 

1.24 

0.147 

1.07 

3.33 

0.019 

1.60 

0.88 

0.39 

1.43 

0.398 

1.63 

0.206 

1.69 

5.21 

0.024 

1.70 

0.88 

0.46 

1.78 

0.510 

2.52 

0.330 

3.26 

1.99 

0.029 

1.77 

0.88 

0.53 

2.14 

0.615 

3.66 

0.476 

5.67 

1.07 


We also note that the effective masses only increase moderately with increasing density. 
However, they are fairly constant over the density range of interest. As a conseqnence, we 
can identify the density n as the major component driving the increases in reTln(2Tp/T) 
and r]T'^ for P = 0 and V = 1. 

On the other hand, for the mixture film in Table HVl DT^ In (2Tp/T) follows an irregular 
pattern with increasing density. This is primarily due to the drastic variation of |A|qP as 
can be seen in column 4 Table |Vl This behavior is due to the fact that is calculated 
from Eq. fl4.4p : = A^q — A^q. Thus a small change in the difference between the two 

forward scattering amplitudes can result in a significant change in the backward scattering 
amplitude. We note that in comparing the denominators of the mixture film transport 
coefficients |Aj -|- |Aj^^Qp -|- |Aj and thus we expect for example that the spin 

diffusion coefficient for the mixture films may be the most sensitive quantity with regard to 
our use of the lowest order ‘T = 0” approximation for numerical calculations. 

In Figs. [5l El and[7]we show the polarization dependence of the thermal conductivity, shear 
viscosity times temperature squared, and the spin diffusion coefficient times temperature 
squared, respectively, for ^He on graphite at re = 0.0252 A“^, and ^He on a 4.33 A ^He 
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TABLE V. ^He in a 3.14 layer film of '^He. The dimensionless scattering amplitudes q’ 
and Landau parameter Fq that are the input into calculating the transport coefficients shown in 
Table HVl 


Density (A 

Obt + iTi.1' + ipot 

V = 0 

FS 

v = o 

II fP 

0 to 

Ipop 

V = 1 

0.013 

0.20 

-0.11 

0.0098 

0.71 

0.016 

0.16 

-0.08 

0.0036 

0.75 

0.019 

0.14 

-0.07 

0.0022 

0.76 

0.024 

0.14 

-0.12 

0.0063 

0.76 

0.029 

0.14 

-0.17 

0.0132 

0.77 


film at n = 0.0248 A“^. The data for k, were calculated at a temperature T = 5 mK which 
was chosen to ensure that the inequality T < Tp 4 , is obeyed at all polarizations. The units 
for K are different in the figure than in the table because for this quantity the temperature 
dependence is not factorable for 0 < P < 1. The major prediction for this section then is 
that K and rj increase dramatically, by an order of magnitude, as V increases from 0 to 1 for 
both substrates. The spin diffusion coefficient goes through a similar large increase from its 
zero-polarization value to its maximum value in the region V ~ 0.74 for both graphite and 
'^He, and then vanishes in the full polarization limit. D vanishes like (1 — P)^P in the limit 
of full polarization (see Sec. IIII Cll . Thus we predict an increase in D from zero polarization 
to its maximum value of 1 ~ 2 orders of magnitude. 

In a recent interesting development, Kovtun, Son, and Starinets^S have conjectured that 
there is a universal lower bound to the ratio of the the shear viscosity to entropy density: 


4^4A>1. 


s/k 


B 


(4.7) 


where s = S/A, and S is the entropy. The authors describe the distance from the lower 
bound as a way to characterize how close a fluid is to being perfect. They argue that possible 
systems that may satisfy the lower bound ought to be strongly interacting systems that are 
normally characterized by a small viscosity (i.e. a small mean free path). They suggested 
that quark-gluon plasmas, and ultra-cold gases at the unitarity limit are candidates. There 
is now evidence that an ultra-cold Fermi gas nearly satisfies the lower bound.— Using the 
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TABLE VI. ^He in a 4.33 layer film of ^He. The zero-polarization 7^ = 0 and full polarization 
V = 1 thermal conductivity k, shear viscosity ry, and spin diffusion coefficient D as functions of 
areal density re with the explicit temperature dependencies factored out. The units of kT are (10“^ 
erg s“^) ; the units of ryT^ are (10“’^ g s“^ mK^); the units of DT^ are (10^ cm^ s“^ mK^). For 
this mixture film niH = 1.29rre is the hydrodynamic effective mass.— The values of re are obtained 


from (I3.45P and (j3.47p . for ry from (I3.80p and (|3.95p . and for D from (I3.120p . for zero and full 
polarization respectively. 


Density (A 

m* 

V = 0 

/mn 

V = 1 

(-F 

V = 0 

(K) 

r = 1 

reT ln(2rp/T) 

r = o V = i 

C-H 

o 

II 

■2 

V = 1 

ln(2rp/T) 

V = 0 

0.015 

1.22 

0.83 

0.50 

1.46 

0.68 

2.29 

0.28 

1.88 

1.74 

0.019 

1.32 

0.84 

0.56 

1.77 

0.80 

3.03 

0.40 

3.03 

2.05 

0.022 

1.37 

0.84 

0.62 

2.03 

0.93 

3.87 

0.53 

4.46 

1.93 

0.025 

1.40 

0.84 

0.70 

2.34 

1.08 

5.18 

0.72 

6.89 

1.61 

0.028 

1.45 

0.84 

0.76 

2.63 

1.21 

6.40 

0.91 

9.60 

1.59 

0.031 

1.50 

0.84 

0.82 

2.92 

2.17 

7.73 

1.08 

12.9 

1.42 


TABLE VII. ^He in a 4.33 layer film of ^He. The dimensionless scattering amplitudes qj 
and Landau parameter Fq that are the input into calculating the transport coefficients shown in 
Table EB 


Density (A 

iTioi^+iVoi^+ 

■p = 0 

II 

14'^'*-12 

V = 0 

iVoP 

P = 1 

0.015 

0.23 

-0.16 

0.03 

0.64 

0.019 

0.21 

-0.16 

0.02 

0.68 

0.022 

0.21 

-0.18 

0.03 

0.71 

0.025 

0.22 

-0.21 

0.04 

0.71 

0.028 

0.22 

-0.22 

0.04 

0.72 

0.031 

0.22 

-0.24 

0.04 

0.73 
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FIG. 5. The thermal conductivity k (13.411) as a function of polarization for ^He on substrates of 
graphite (solid line), and a 4.33 A superfluid ^He film (dashed line). Both results are at T = 5 mK, 
and the ^He areal densities are 0.0252 and 0.0248 A“^ on graphite and ^He, respectively. 


results from Sec. IIII Bl we can estimate the value of this expression for a ^He him. From 
Fig. [6] we see that the viscosity is a monotonically increasing function of polarization. The 
polarization dependent entropy density is given by s/ks = (vr/Oh^) (m| + m|) ksT.— This 
entropy is a monotonically decreasing function of polarization. Thus, we need only to concern 
ourselves with the zero-polarization limit. Using fl3.80p we hnd for the left hand side: 


dvr 


dvr 


s/ks 

v/^ 

s/k 


(m/m*) 


* \ 2 


TT^ 


I |2 I A^l 12 I dtt _ dtt 12 


U,0 


ffi 


Uoi 


T 


B 


0.28 l'^ 


(4,8) 


where we have used (ld.2p to write the transition rates in terms of the dimensionless scattering 
amplitudes. The numbers come from Table [TTl and Fig. 01 and so they refer to the second 
layer of ^He on graphite at rT = 0.025 A~^ (Tp = 0.7d K). It is clear from the inverse cubic 
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FIG. 6. The shear viscosity r] (13.941) times temperature squared as a function of polarization for ^He 
on substrates of graphite (solid line), and a 4.33 A superfluid ^He film (dashed line). The results 
are shown for ^He areal densities of 0.0252 and 0.0248 A“^ on graphite and ^He, respectively. 


temperature dependence that deep in the Fermi-liqnid regime the system satishes the lower 
bonnd. At higher temperatnres this expression passes throngh one when T 0.5 K. This 
is not that high, and snggests that at temperatures on the order of lOO’s mK the ratio may 
not be very far from one for this ^He thin him system. 


V. CONCLUSION 

We have derived exact expressions for the transport coefficients k and r] utilizing methods 
developed by numerous groups^ for application to bulk ^He. We calculated predicted values 
for the polarization dependence of k, t], and D for thin, degenerate ^He hlms using previously 
determined Landau parameters. The key to performing the principal angular integration 
in phase space is the procedure suggested by Miyake and Mullinl^ for avoiding a hnite 
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FIG. 7. The spin diffusion coefficient D (|3.114p times temperature squared as a function of po¬ 
larization for ^He on substrates of graphite (solid line), and a 4.33 A superffuid ^He film (dashed 
line). The results are shown for ^He areal densities of 0.0252 and 0.0248 A“^ on graphite and 
■^He, respectively. 


temperature singularity. The Miyake-Mullin approach is discussed in detail in Sec. HIl In 
that section we derive the polarization dependent expression for the quasiparticle lifetime due 
to quasiparticle-quasiparticle collisions in the relaxation time approximation. We compare 
that result with that of a previous derivation of the quasiparticle lifetime using completely 
different techniques, and note that they are identical up to factor of order one. 

The derivation of the transport coefficients in Sec. IIIII follows the methods developed by 
Abrikosov and Khalatnikov,- and Sykes and Brooker.— The calculation of k is very similar to 
that of the spin diffusion coefficient D as described by Miyake and Mullin. The collision in¬ 
tegral is reduced to an integral eigenvalue problem whose integrand depends on both spin-up 
and spin-down fluctuations. The system is diagonalized by standard methods, and is reduced 
to an independent pair of equations in Sykes-Brooker form. The temperature dependencies 
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TABLE VIII. The temperature dependencies of the inverse transport coefficients as a function of 


polarization. The thermal conductivity k 
spin diffusion coefficient D is from Ref. [h 

and the shear viscosity rj are calculated in 
]. We note that D~^ is undefined at R = 1, 

Sec. Uni The 

Coefficient 

o 

II 

0 < R < I 


TlnT 

TlnT 

jy-i 

T2 

T2 


InT 

rjp2 


for the transport coefficients are in agreement with older work at zero polarization by Fu and 
Ebner.— Further, we find that, unlike spin diffusion, these dependencies (TlnT for k and 
for rj) are not changed by polarization. The solution for the shear viscosity is unlike that of 
any other fermion transport coefficient. The key physics lies in including the contributions 
of scattering from quasiparticles whose momenta differ slightly from their zero-temperature 
values but are still allowed by energy and momentum conservation at non-zero temperature. 
We introduced a simplified model in which we fix the incoming quasiparticle momenta at 
the zero-temperature values, and allow the outgoing momenta to vary (see Fig. [3]). We find 
that in lowest order the viscosity is formally independent of the quasiparticle lifetime (see 
Eq. (13.791 for example). We note however that I/z/q (13.781) is very similar to Tq. In Ref. 


24 


Novikov, in the zero-polarization limit, allows all four quasiparticle momenta to drift from 
their zero-temperature values. We find at zero polarization, in agreement with Novikov, that 
the head-on collisions between quasiparticles with momenta in opposite directions dominate 
the scattering process, and we also find that they are the dominant process in the scattering 
between spin-parallel quasiparticles at finite polarization. Our final result for the shear vis¬ 
cosity temperature dependence differs from that of Novikov because Novikov assumes the 
Landau parameters have a divergence at 6* = tt, and this gives an extra factor of In^ (Tp/T) 
in the final result for the viscosity. 


In Sec. IIVI we apply these results to a system of thin ^He films both in the second layer 
on a graphite substrate, and also in a thin ^He-'^He film mixture. In Table IVIIII we gather 
together the main results from this paper concerning the temperature dependence of the 
thermal conductivity and the shear viscosity, and we have also included the spin diffusion 
coefficient results from Miyake and Mullin.— 
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The predicted polarization dependence of the transport coefficients for ^He on the second 
layer of graphite and also for the 4.33 A-thick ^He film is shown in Figs. O [6l and [71 These 
results show a dramatic increase in the magnitudes of the coefficients as the polarization 
increases from zero. We showed in Sec. IIII Al that for the thermal conductivity in two 
dimensions k is proportional to the quasiparticle lifetime. Further, we showed in previous 


work, see Fig. 7 in Ref. [l^ that the magnitude of the contribution to the quasi-particle 


lifetime from the majority spin component decreases dramatically as a function of increasing 
polarization. Thus, for the thermal conductivity a fairly simple qualitative picture emerges 
of the role of polarization: increasing V induces an increase in the quasiparticle lifetime, 
and thus the transport coefficient. For very dilute systems this mechanism is basically 
understood as the quenching of s-wave scattering with increased V. However for the shear 
viscosity such a simple picture does not seem to be relevant if for no other reason than 
because the quasi-particle lifetime does not contribute directly to the transport coefficient. 
In this case we must consider instead the complicated dynamical question of the relative 
importance of the spin anti-parallel head-on scattering to the spin parallel head-on scattering 
as per the discussion in Sec. IIII B1 which itself is related to the balance of s-wave and p-wave 
scattering. 

In lowest order of temperature the derivation of the expressions for the transport co¬ 
efficients is essentially exact. The calculation of explicit results for ^He films suffers from 
the use of the £ = 0 approximation for the scattering amplitudes. An improvement in the 
present results would be the inclusion of additional Fourier components in the expressions 
for the transition probabilities in terms of the scattering amplitudes, see Eq. fl2.22p . The 
approximations used in the determination of the ^He film Landau parameters from exper¬ 
imental measurements of the specific heat effective mass, and the spin susceptibility have 


been discussed in Ref. 
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At this time to the best of our knowledge there have been no measurements of any 
transport coefficient in a thin ^He film. In addition there have been no measurements at 
all in a polarized thin ^He film. These experiments would be very difficult. In fact the 
first measurement of zero sound in a thin, unpolarized ^He film was only reported in 2010 
by Godfrin, Meschke, Lauter, Bohm, Krotscheck, and Panholzer.— Our Landau parameters 
do yield excellent agreement with this zero sound measurement. For bulk ^He there has 
been some work on the polarization dependence of transport coefficients. A recent review^ 
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summarizes the state of the held. 
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